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A semi-analytical study of the acoustic radiation losses associated with various transverse vibration 
modes of a micromechanical (MEMS) annular resonator is presented. The quality factor, Q, of 
such resonators is of interest in many applications and depends on structural geometry, interaction 
with the external environment, and the encapsulation method. Resonators with at least one surface 
exposed to air can display losses through acoustic radiation even at fj,m dimensions. Published 
experimental results suggest the dominance of acoustic losses in the Q of a MEMS drum resonator. 
In this study, a well established mathematical techniques to analytically model resonator vibration 
modes and fluid-structure interaction are used, and a semi-analytical procedure for computing Q 
due to acoustic radiation losses, Qac, in any vibrational mode outlined. Present technique includes 
calculation of the exact mode shape and its utilization in computing Qac. The dependence of Qac on 
the first 15 mode shapes is computed. Results are compared for the lowest 2 modes of a solid circular 
resonator using exact mode shapes to those of Lamb's approximate mode shapes. Comparison to 
published experimental results validates the predictive utility of the technique, especially for higher 
modes where acoustic radiation seems to be the dominant constituent of Q. 



PACS numbers: 4340Rj, 4340Dx, 4340Ey, 4338Gy 
I. INTRODUCTION: 

Mechanical resonators can have significant advantages 
such as small size, low power operation^ and the abil- 
ity to be integrated into CMOS structures^ yielding high 
quality factor (Q) devices 3 for use as RF filters 4 and for 
mass sensing^. For chemical sensing, the route that uses 
changes in stress to monitor the presence of trace vapors 
from the ambient air is particularly attractive^^. In 
this sensing method, a polymer coating applied to one 
side of the mechanical resonator changes its morphol- 
ogy, swelling or contracting and thus altering the stress 
in the underlying silicon structure, which in turn pro- 
duces a shift in the resonant frequency. However, the 
fact that the vibrating mechanical structure must be im- 
mersed in air so as to enable the sensing mechanism has 
consequences, since the presence of the air leads to me- 
chanical energy loss to the medium. The loss increased 
dissipation Q _1 , is manifested as a broadening of the 
resonant peak, leading to a degradation of the ability to 
resolve frequency shift and consequently reduction in the 
sensitivity to analyte concentrations. 



a 'Electronic address: santdmec@mecheng.iisc.crnet.in Also at De- 
partment of Mechanical Engineering, Indian Institute of Science, 
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It has long been recognized that there are two main 
sources of energy loss of such micron scale resonators. 
The first, termed squeeze film losses are due to the 
trapped film of air under a vibrating resonator^ 1 ^. The 
second loss mechanism is where the vibrating structure 
couples acoustically to the air. This paper quantifies the 
loss due to acoustic radiation for a conceptually simple 
readily fabricated structure: an annular plate fixed at its 
edges. We show that the acoustic losses can be reduced 
by operating the structure at higher harmonics, where 
the phase difference between adjacent moving segments 
at various vibrating modes of the resonator leads to a 
reduction in radiated energy. 

A circular plate vibrating in contact with surrounding 
fluid, is of significant interest in wide range of systems 
such as the piston movement in contact with fluid in a 
closed cylinder—, ocean structure moving underwater—, 
movement of magnetic disk drive in contact with air—, 
vibration of nano circular drum resonator in contact with 
surrounding fluid in case of pressure sensors 1 ^, microflu- 
idic device with circular membrane^, etc. There are 
two main effects associated with such systems due to the 
coupled effect of fluid-structure interaction. In the first, 
the "added mass effect" reduces the effective resonance 
frequency. In the second, the acoustic radiation affects 
the quality factor of the vibrating structure. This fluid- 
structure problem is analyzed numerically using the finite 
element method or boundary element method for com- 
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plex domains and analytically for simple domains. In this 
paper, we present an analytical procedure to analyze the 
acoustic radiation losses in higher modes of the vibrating 
annular plate with fixed outer boundaries. 

Lord Rayleighi& was probably the first to study ana- 
lytically the effect of the mass loading of the surround- 
ing fluid on a vibrating rigid disk in contact with the 
fluid and suggest the idea of "added mass". This clas- 
sic problem was then studied by Lamb 12 to investigate 
the added mass induced change in the first two natural 
frequencies of a circular plate, fixed along its outer edge 
and vibrating in contact with a fluid. In addition to the 
frequency change, he also found the acoustic radiation 
losses corresponding to those two modes. These results 
have been validated experimentally by subsequent stud- 
ies. It is, however, important to note that Lamb's results 
are based on approximate mode shapes. Revisiting this 
problem, Amabili et a L 17 ' 18 used the Hankel transform 
to analyze the added mass effect on the frequency of the 
structure vibrating in its fundamental and higher modes. 
However, their study was limited to analysis of the added 
mass effect. Such fluid-structure interaction effects have 
recently been analyzed in the context of Micro and Nano 
Electro Mechanical Systems (MEMS/NEMS) structures. 
Unlike in the case of macroscale problems where the sur- 
rounding fluid is assumed to be incompressible and in- 
viscid, the experimental validation of Lamb's theory was 
done examined for a microsensor working in the pres- 
ence of a viscous fluid (water-glycerol mixture)^. Ex- 
perimental results were found to be in good agreement 
with Lamb's predictions for a less viscous fluid mixture 
(< 10 cP) but differed for fluid mixtures with higher vis- 
cosity. The difference was attributed to the viscosity con- 
tribution to the added mass effect qualitatively^ as well 
as quantitatively 2 ^. Considering the first approximate 
mode shape of a circular plate vibrating in contact with 
the surrounding fluid, Kozlovsky 20 analyzed the effect of 
viscosity on the natural frequency as well as the quality 
factor. Recently, Olfatnia et al^ have compared theo- 
retical and experimental results for analyzing the effect of 
viscosity of the surrounding fluid on the frequency as well 
as the quality factor of a circular diaphragm vibrating in 
its first mode. 

Our central interest in the present study is to find 
acoustic losses and the associated quality factor in var- 
ious modes of vibration of an annular micromechanical 
resonator (a MEMS plate) clamped at its outer edge in 
order to assess the suitability of higher modes to high- 
Q applications. The annularity of the resonator results 
from the requirement of an etch hole typically used in 
the micromachining technique to create a cavity under- 
neath the resonator by etching away the supporting ox- 
ide material starting from the etch hole. Although the 
etch holes can be many and spread over the plate, a cen- 
tral etch hole serves the purpose and creates the simplest 
resonator structure^. The annular resonator is also par- 
ticularly amenable to analytical treatment for studying 
acoustic losses and building mathematical models for pre- 
dicting the Q-factor. 

In order to analyze acoustic radiation losses associated 
with various modes of vibration of the annular plate, we 



first derive the exact mode shapes of the structure ignor- 
ing any effect of the surrounding fluid (air) on the mode 
shape. We use these mode shapes to study the effect 
of the surrounding fluid on the associated natural fre- 
quencies and the Q-factor. Since the surrounding fluid 
is air, the effect of "added mass" on the frequencies of 
the structure is negligible. The effect of the surrounding 
air on the Q-factor, however, is significant because of the 
acoustic radiation losses which is the subject matter of 
this paper. We extend the analytical approach proposed 
by Amabili et aL 17 ' 18 to higher modes, establish an ana- 
lytical calculation procedure involving symbolic algebra, 
find the acoustic radiation losses, and compare our re- 
sults, first, with Lamb's results for the first two modes, 
and then some published experimental results^ for the 
higher modes of the resonator. 



II. MATHEMATICAL MODELING 

In this section, we present the mathematical back- 
ground for computing the acoustic radiation losses due 
to a vibrating annular plate of free inner edge and fixed 
outer edge. The annular plate, with the inner radius a 
and the outer radius b in contact with the hemispherical 
surrounding fluid on its upper surface as shown in Fig. [T] 
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FIG. 1. Schematic of semi-hemispherical fluid domain above 
the vibrating annular plate fixed at its outer edge 



First, we find the expression for the exact mode shapes 
of the vibrating annular plate in vacuum. Subsequently, 
using the same mode shapes, we present the procedure 
for computing the acoustic radiation losses. 



Acoustic Radiation Losses of Micro Resonators 2 



A. Modal Analysis of Annular Plate 

The vibration of annular plates with fixed outer bound- 
ary is widely analyzed for different applications. Under 
the assumptions of a thin plate made of isotropic, ho- 
mogeneous, and linearly elastic material, the equation 
governing the transverse deflection W(r, 9, t) of the plate 
for undamped, free vibration can be written in polar co- 
ordinates as ) 22 ' 23 i 24 



V 4 IU - 



Ps d \ (PW_ 
dt 2 



D 



0, 



(1) 



where, the Laplacian operator V 2 = + --^ + -^Jjt 
and D = ln ^ d °..i\ is flexural rigidity of the plate, p s is 



12(1-1^) 



the density and do is the uniform thickness of the plate. 
Assuming that the plate vibrates in a normal mode, W 
can be expressed as, 



W = u;(r,0)e* 



(2) 



Substituting eqn. ^ in eqn. (fTJ) and then rearranging 
the resulting equation, we get the modal equation as, 



V 4 w(r, 9) - (3 4 w(r, 9) = 0, 



(3) 



where, /3 4 = ^nw 2 , y is the frequency parameter. The 
modal solution w(r, 9) can be expressed in terms of Bcsscl 
functions as described in Appendix (|A"]) and is given by^SS 



system of linear homogenous equations (|B1[) - (|B4[) . Since 
the equations are homogenous, these constants can be ex- 
pressed in terms of any one (say D mn ) of these constants. 



Therefore, expressing the constants A 



null B mn . 



terms of D mn and substituting them in eqn. (|A10|) . the 

mode shape can be written as ^fey^ = D mn T mn ((3r). 

Normalizing the mode shape based on the normalization 
condition^ 



j(a,9) 



ada = 1, 



ib/a \ m 

where, a — r/a, we get the expression of D mn as 



\J fl{T mn (f3aa)} 2 ada 



(5) 



(6) 



For a particular resonant mode of vibration of the annu- 
lar plate, the computed D mn from eqn. (J6j) can be used 
to determine other constants A mni B mn and C mn respec- 
tively. 

For the limiting condition of the annular plate when 
the radius of the inner hole goes to zero, i.e., a solid 
plate, we find, following a similar procedure as described 
above, that Y n and K n tend to infinity as r tends to 
zero 2 ^. As w is finite at the center of the plate, we must 
set B mn and D mn to be zero. The resulting deflection 
profile takes the form ; 24 ! 25 



w(r,9) = x(rM6) = [A mn J n {pr) + B mn Y n (j3r)+ 

C m nln(pr) + D mn K n (/3r)]iP(6), (4) 

where, J n is the Bessel function of the first kind, Y n is 
the Bessel function of the second kind, /„ and K n are the 
modified Bessel's functions of the first and second kind, 
respectively. 

For the annular plate with the fixed outer edge 
and free inner edge as shown in Fig. [TJ the fol- 
lowing boundary conditions can be used. At the 
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ing these boundary conditions, we get the system of four 
linear and homogenous equations for the four constants 



w(r,6) = [A mn J n {/3r) + C mn I n {f3r)]^{9). 



(7) 



Using the boundary conditions and following a similar 
procedure as described above, we get the following char- 
acteristic equation, 



J n (fia)I n+1 (fia) + I n {pa)J n+1 {fia) = 0. 



(8) 



Again, solving the above characteristic equation for /3 
and using the normalization condition, we can obtain ex- 
pression for values of A mn and C mn . 

Most of the fluid structure interaction problems 12 use 
a polynomial approximation for the mode shapes. How- 
ever, such approximation introduces errors in estimating 
the damping from the fluid-structure interaction. In the 
present analysis, we use exact mode shapes for estimating 
the acoustic radiation losses in different mode shapes. 



C mn and D mn which we list as eqns. (|Bl[) - B. Estimation of Acoustic Radiation Losses 



A R 

(|B4p in Appendix [B] For a non-trivial solution of these 
constants, the determinant of the coefficient matrix given 
by eqn. flDlO. Appendix D, of the above equations is set 
to zero, which gives the required characteristic equation 
governing the frequency constant ft. Since the character- 
istics equation given by eqn. ([DID, Appendix D, is dif- 
ficult to solve analytically, we employ a numerical tech- 
nique based on the bisection method for finding roots 
in MATLAB to estimate f3. The different solutions of 
(3 are the natural frequencies, /dry, for different modes. 
For a given value of /3, the corresponding constants A mn , 
Bmn, C mn , and D mn can be determined by solving the 



In the previous section, we found the natural frequen- 
cies and corresponding mode shapes of the annular plate 
vibrating in vacuum. When the plate vibrates in contact 
with the relatively denser surrounding fluid than air, the 
change in the mode shape of the plate is assumed to be 
negligible but there is a non-negligible reduction in the 
modal frequency. Such hypothesis has been used in many 
fluid-structure interaction problem s 17 ' 18 i 27 . However, we 
consider the corresponding changes in frequency in our 
formulation for estimating the acoustic radiation losses 
in the surrounding fluid. 
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As the annular plate vibrates, a disturbance is cre- 
ated in the fluid adjacent to the plate, which causes the 
wave motion to introduce the pressure fluctuations at all 
points in the fluid domain 28 . Considering the surround- 
ing fluid as irrotational and inviscid at constant ambient 
mean pressure po, temperature To and density pf, the 
pressure fluctuation above po and the particle velocity 
can be found in terms of the velocity potential <j) through 
the equation ; 11 1 28 



P 



at 



and v = V0. 



(9) 



The governing equation for the velocity potential corre- 



sponding to small disturbances is given by; 11 ! 28 ' 29 



V 2 



c 2 di 2 ' 



(10) 



where c s is the speed of acoustic waves in the fluid. The 
incompressibility condition V.w = leads to the Laplace 
equation in terms of </>, i.e., V 2 = 0. Therefore, it is 
sufficient to find the velocity potential field in order to 
analyze the propagation of waves in the fluid medium. 

For the domain shown in Fig. [TJ the fixed outer sup- 
port (i.e., for r > a) is assumed to be radially extended to 
infinity. For the hemispherical fluid domain enclosing the 
upper surface of the annular plate and the support, along 
with the Sommerficld boundary conditions, i.e., velocity 
and its gradient vanish as r — > oo, the velocity potential 
can be obtained by solving the wave equation as men- 
tioned above. 

To find the velocity potential at a generic 
point P due to an elementary source at S as 
shown in Fig. [21 we follow the analysis given by 
Lambi 2 -. Let the position co-ordinates of point P be 
(i? sin £ cos sin £ sin cos £), and that of S on 
elemental surface area dS in the plane of the resonator 
be S(r cos 9, r sin 0, 0). The distance r' of the the point 
P from the source position S is given by 



{R 2 - 2 Rr sin £cos(V> ~9)+r 2 }~- 



(11) 



Since at the far field point P, R r, the expression for 
r' can be approximated as, 



r' = {R — rsin£cos(V> — 9)} . 



(12) 



Assuming that the normal component of the velocity of 
the fluid at the plate interface is same as that of the plate, 
the fluid velocity v n = at point S on the upper sur- 
face of the annular plate can be written in terms of the 
mode shape w(r, 9) defined in section Hl.AI When the en- 
ergy loss in the form of acoustic radiation is accounted for 
(damped structural response) , we need to modify eqn. ([2]) 
to accomodate a time varying amplitude that represents 
the transient decay. Assuming W = w(r,9)e lult , we can 



find the normal velocity v n 



i)W 
dt 



as. 



dW 

— — — = tut w(r, I 
dt v 



(13) 



But, w(r,9) = xMVK^) = x( r )C cosn9, (see Appendix 
(A)). Therefore, v n — Ciuix(r) cosn9e lut . 




Clamped Edge 



FIG. 2. Generic point P in acoustic far field from the res- 
onator 



Let v n = Ax(r) cos n9e iUlt , where A = iCco. We also 
know that because of acoustic radiation losses, the am- 
plitude A has to be time dependent, but its time depen- 
dence remains to be determined. So, at this point, we 
assume that, 



v n = A(t)x(r) cosn6 



(14) 



where x(r) = [A mn J„(/3r) + B mn Y n (f3r) + C mn I n (f3r) + 
DjnnKn^r)], uj is the frequency of the vibrating plate 
when in contact with the fluid. The velocity potential 
at a distance R from the center of the resonator, or at 
a distance r' from the elementary surface dS due to the 
disturbance at S is given by ; 12 ' 30 



27r JO 



-ikr 



dW 



dS. 



(15) 



Using the approximation of r' from eqn. (|12p and the 
expression of ^f- from eqn. ([M)) . the velocity potential 
can be rewritten as, 



A(t) 



,ik{c t t-R) 



2nR 



271 



^ikrsin£ cos(t/j — G) 



x cosrt6'x(r)rdrd6' (16) 



where k = — is the wave number. Using the following 
property of Bessel functions , 13 i 29 

r 2ir 



^ikrsin^ cos(Tp — 9) 



cosn9d9 = 2ni n cos nipJ n (kr sin £), 



o 



the velocity potential becomes, 
A ^ ine ik(c,t-R) 



(17) 



R 



cosmp x / J n (kr sin t;)x(r)rdr. 

(18) 
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From this expression of velocity potential, one can calcu- 
late the pressure fluctuation and velocity using eqn. ([9]). 
Taking the real part of the velocity potential, 



= Re 



(A{t) 



;n ik(c s t — R) 



R 



cos nip 



x / J n (kr sin £)x(r)rdr, (19) 

Jb 

corresponding to the real part of the velocity ^f- = 
x(r) cos (n0) cos (kc s t), the intensity of the acoustic wave 
at any point at distance r' can be written as the product 
of pressure and particle velocity at that poinljii. The cor- 
responding power radiated across a hemispherical surface 
of radius R, is obtained by integrating the intensity over 
the hemispherical fluid domain as 



dE Ru 
dt 



271 



JO 



-PS 



^i^W (20) 



dt OR 



The mean flux of energy emitted in the form of sound 
waves is given by 



E 



f W WCt dEflvLX (If 

aii? = ^_ dt - = UA 2 (t). 
/ — dt 



(21) 



The expression for U is obtained by substituting the ve- 
locity potential, ip r in the power radiated eqn. (|20[) and 
performing time averaging as in eqn. (|2ip . The parame- 
ter U depends on fluid properties, geometric properties 
and frequency of oscillation of the resonator and is given 
by, 



U = -p f R- 



27r Jo \Jb 



J„(fcrsin£)x(r)rdrJ sin£d£ 
x / cos 2 (™/>)oV I "' -^(Re(z n e lfe(Cst - it) ) 



2- g 



dR V R 



dt. (22) 



The total energy pumped into the system by the source 
is the mean kinetic energy of the plate and the adjacent 
fluids 



£input =T P (1 + Pmn) = VA 2 (t), 



(23) 



where Tp is the mean kinetic energy of the resonator, 
and (3 mn is the added virtual mass incremental factor- 
as mentioned in appendix [E] The mean kinetic energy 
Tp is given by, 



Tp = ^p s d 



2tt pa 



Jb 



(x(r)A(t) cos n9) 2 rdrd9 



V 



2(l + /3 mn ) 



A 2 (t). (24) 



On equating the rate of decay of input energy from the 
source, i.e., 



dE iuput _ y A{t) dA (t) 



dt 



dt 



(25) 



with the rate of energy emitted into acoustic field, i.e., 
eqn. (|2"T]) . we gei^. 



dA{t) 
dt 



= SA(t), 



(26) 



where <5 = U/V, the closed form expression for S is given 
by, 

pj R f wet Jo 5 (J b a J„(fcrsin^)x(r)rxir) 2 sin^dC 



p s d (1 + f3 mn ) 



(c s t-fl)- 



J b a X 2 (r)rdr 

d /Re(i ^^ e ^fc ( c = t -• R) ) , 
dR\ R 



(27) 



Equation (1261) shows that the amplitude diminishes ex- 
ponentially with an exponential constant 5. On equating 
this constant 5 with the decay constant, C^wct of the 
damped linear oscillator—, we get C = S/uj wct . Based on 
the definition of quality factor—, we get 



Qa 



1 

2( 



T/dry 



(28) 



where u; wc t — 27r/ wct , / we t is the frequency of the vibrat- 
ing plate in contact with the fluid. Although the final 
expression of the quality factor looks very simple, it re- 
quires the computations of 5 and j3 mn which are mode 
dependent. Since the computation of these parameters 
requires successive differentiation and numerical integra- 
tion, we perform these steps in MATLAB and MAPLE. 
The entire set of calculations is somewhat involved but 
algorithmic. We, therefore, present the sequence of cal- 
culations to be done as a schematic flow chart in Fig. [3] 
There are three major computational blocks that we 
have tried to present separately for conceptual clarity. 
As is evident from Fig. [3l the computation involves one 
FEM analysis along with computer algebra and general 
numerics. 

We can also use the procedure described above to com- 
pute the quality factor for a solid circular plate by setting 
6 = and x(r) = A mn J n (/3r) + C mn I n (/3r). In the sub- 
sequent section, we compare the quality factor obtained 
using exact mode shapes as described in this paper with 
that using approximate mode shapes based on Lamb's 
formulation^^ for the first two modes of plate. Next, we 
use the analysis presented here to compute the quality 
factor for the annular resonator in higher modes, covering 
15 modes. We also compare the computed quality fac- 
tor with experimentally obtained values^ in the higher 
modes where acoustic radiation dominates the losses. 



III. RESULTS AND DISCUSSION 

We first compute the quality factor for a clamped solid 
circular plate vibrating in contact with the surrounding 
medium. We concentrate on the first two modes. We 
compare our results with those of Lamb. 

To do the analysis, we take the following dimensions 
and material properties of the structural and the fluid 
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Common Input Block 
Structural geometry and material properties 
a, b, d , q,E, v, fj, c s 
(used in all computational blocks) 



kg/m 3 under constant ambient temperature T = 293 K 
and pressure P = 1.013 x 10 5 Pa. For the analysis of the 
solid plate, we simply take the inner radius to be zero, 
i.e., b — and rest of the parameters remain the same. 



Modal Analysis 

(following Vogel and Skinner [22]) 

• Carry out modal analysis (ANSYS) to get all frequencies 
in the range of interest 

• Use f di corresponding to the mode '« ' of interest as the 
initial guess and compute initial value of (3 used in eqn.(3) 

• Solve for exact [3 from eqn. (Dl) using numerical 
iteration (e.g., MATLAB calculation using bisection 
method) 

• Using J3 and 'h ', compute exact f, and the correspond- 
ing mode shape co-efficients Arm, Bum, Cum and Dmn, 
from eqns.(Bl-B4), (5) and (6) 



Computational 
Block - : 



INPUT 



Fluid Structure Interaction 

(following Amabili [17,18]) 

• Formulate functions in eqns. (E4-E12) 
(MAPLE) 

• Using appropriate limits of the integral in 
eqn. (E3), split the integral and evaluate 
using quadrature ( we use Gauss-Kronrod 
quadrature ) to obtain [3 



Computational 
Block - 2 



OUTPUT 



Acoustic Radiation Loss 

(following Lamb [12]) 

• Formulate velocity potential, (J) given in 
eqn. (18) (MAPLE) 

• Substitute potential in eqn. (20) and find the 
mean flux of energy E from eqn. (21) and get 

MF 

an expression lor U 

• Compute mean kinetic energy of the plate, T p 



expression for V. 



• Compute § = 



Computational 
Block - 3 



, OUTPUT 



INPUT 



1 



Obtain 0* 



FIG. 3. Flow chart for computing acoustic radiation damping 



A. Solid circular plate 

We compare the frequency and quality factor based on 
the exact mode shapes given by eqn. (J7J) and the approxi- 
mate mode shapes used by Lamb for the first two modes. 
Fig.[4]and[5]show the comparison between the normalized 
exact mode shape (ems) and the normalized approximate 
mode shape (ams) for the axisymmetric mode shape with 
zero nodal circle, i.e., (to = 0,rt = 0) and single nodal 
diameter mode (to = 0, n = 1), respectively. As is evi- 



Exact (e) Circular Mode 
Approximate (a) Circular Mode 




hi 
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FIG. 4. (color online) Normalized circular mode shapes of a 
clamped solid circular plate 



dent from these figures, the exact mode shape differs from 
the approximate one in spatial amplitude variation. The 
variation is significant enough to cause non-negligible dif- 
ferences in the radiated acoustic power. This is precisely 
what we see in the computed values of corresponding Q, 
listed in Table HI The difference in amplitude variation 
across the plate between the two mode shapes causes 
very significant change in the corresponding Q. In par- 
ticular, the approximate mode shape used by Lamb in 
the second mode overestimates the radiative losses by as 
much as 42%. This approximation may lead to even big- 
ger difference when we try to compute the Q's for an 
annular plate in the next section. 



domains based on the test structure used by Southworth 
et al*£ in their experimental studies. Here, the annular 
plate is made of silicon with Young's modulus E = 150 
GPa, the density p s = 2330 kg/m 3 , and the Poisson's 
ratio v = 0.22. It has an outer radius of a — 18.4 /xm, 
inner radius of b — 2 /im and a thickness of do = 300 nm. 
The plate is surrounded by air having density pj = 1.2 



B. Annular plate 

We now consider the annular plate shown in Fig.[T] We 
intend to compare our analytically computed results with 
experimental results of Southworth et a/J^. Therefore, 
we take the geometry and the material properties used 
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mode of vibration. Fig. [6] shows the computed values of 
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o 

o 
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o 




• • 
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o 

Acoustic radiation 


o 

o 


dominated region 

o 



Exact (e) Diametral Mode 




Normalized Radius (r/a) 

FIG. 5. (color online) Normalized diametral mode shapes of 
a clamped solid circular plate 



TABLE I. Quality factors associated with exact and approx- 
imate mode shapes corresponding to circular and diametral 
modes 



m 


n 


/wet 


Qcms 


Qams 








(MHz) 














3.38 


116 


100 


13 





1 


7.08 


233 


135 


42 



in Southworth et al. for their drum resonator. For this 
resonator, the annular plate has hxed outer edge with the 
ratio of inner to outer radii as 0.1. Moreover, the lower 
surface of the plate is close to another fixed-fixed plate 
but the upper surface is open to the surrounding air. 
Based on this geometric configuration, there are two ma- 
jor sources of energy dissipation: the squeeze film damp- 
ing due to the trapped air film in the cavity below the 
resonator and the acoustic damping due to the acoustic 
radiation in the free space above the resonator surface. 
We use the procedure outlined in section II to obtain the 
exact mode shapes of the resonator and use these mode 
shapes to compute the corresponding Q for each mode. 
The results are shown in Fig. [TTJ We also include the ex- 
perimental results reported by Southworth et al. in this 
table. Our focus here is on computing Q ac - the qual- 
ity factor associated with acoustic radiation losses - and 
understanding its variation in different modes of vibra- 
tion. Because acoustic radiation is not the only damping 
mechanism, we do not expect to match the experimen- 
tally obtained Q values with the computed Q ac values 
here. A comparison of the two Q's listed in Fig. HT1 clearly 
indicates that the relative contribution of Q ac is very low 
in the lower modes of vibration and becomes dominant 
in the higher modes. It is well known that squeeze film 
damping is high at lower frequencies^ and hence proba- 
bly contributes maximum to the damping up to the 4 th 



2 4 6 8 10 12 14 
Mode No. 

FIG. 6. (color online) Computed acoustic quality factors for 
various modes 

Q ac graphically. It is interesting to see that the variation 
in Q ac is not monotonic with increasing frequency over 
the entire frequency range. It has a markedly different 
behavior in the lower modes (upto 5 th mode), but settles 
down to a more predictable, almost monotonic trend, in 
the higher modes. The unsettled behavior in the lower 
modes can be explained qualitatively by looking at the 
corresponding mode shapes in Fig. Inland pondering over 
the relative efficiency of these modes in transferring the 
resonator's energy into the acoustic radiative field. As 
the mode number increases, the resonator starts behav- 
ing like many point sources and addition of more such 
sources with increasing mode number makes less and 
less difference to Q ac leading to a somewhat predictable, 
monotonic and slower increase in Q ac - 

We know that l/Q no t = 1/Qac + 1/Qrest- We use 
this relationship to evaluate the relative contribution of 
Qac to Q nct (which is the Q oxp here). The last col- 
umn of Fig. [H] lists this contribution as a percentage of 
Qnet (computed from Q GX p/Qac because of reciprocal re- 
altionship). This column clearly shows the dominance of 
Qac in higher modes. Fig. [7] is a graphical representa- 
tion of Q ac and Q cxp in the acoustic radiation dominated 
modes. This figure shows that Q ac and Q oxp are close 
to each other in these modes and thus the net Q is pre- 
dominantly dependent on the acoustic radiation. Hence 
any effort to enhance Q in these modes should focus on 
decreasing acoustic radiation losses. 

We also investigate the influence of the ratio of in- 
ner and outer radii on the quality factor. We observe a 
marginal increase in the acoustic quality factor with in- 
creasing radii ratio for all representative resonant modes 
shown in Table ILU1 It is not until b/a = 0.7 (representing 
about 1/2 of the plate area as the hole) that we see a 
significant change in Q ac , particularly in lower resonant 
modes (see Fig. [5]). At higher modes, the increase in Q ac 
is marginal even for this large b/a ratio. This is along ex- 
pected lines as the resonator radiates less and less acous- 
tic power in the higher modes and any reduction in the 
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TABLE II. Experimental and theoretical comparison of acoustic Q-factors for annular resonator clamped at the outer edge(a) 
and free at the inner edge(fr) with b/a = 0.1087 



Mode No. 




Y n (Pr) + Cmnln(Pr) + D mn K n (Pr) 


AVMI 


fwet 


fexp 


Quality Factor 


% Share 


(m, n) 


P 




Bfnn 


Cm n 


Dmn. 


{3mn 


(MHz) 


(MHz) 


Qac 


Qcxp 


Of Qac 


1(0,1) 


2.50149e5 


3.84177 


-0.03135 


0.05894 


0.03846 


0.09407 


7.06 


5.53 


236.42 


16.84 


7 


2(0,2) 


3.18878e5 


-4.43796 


0.15639 


-0.02419 


0.10691 


0.00656 


11.49 


9.79 


260.02 


44.22 


17 


3(0,3) 


3.87993e5 


-5.02600 


0.01979 


-0.01016 


0.01377 


0.00515 


17.02 


14.95 


206.13 


74.61 


36 


4(1,1) 


4.20944e5 


4.89110 


-0.26138 


-0.00314 


0.40224 


0.00334 


20.05 


21.08 


214.78 


95.38 


44 


5(1,2) 


4.95634e5 


-5.37671 


0.40439 


0.00102 


0.27212 


0.00415 


27.79 


28.05 


161.15 


138.32 


85 


6(1,3) 


5.71687e5 


-5.87366 


0.09854 


3.38038e-4 


0.07402 


0.00346 


36.98 


32.97 


208.94 


176.67 


85 


7(0,7) 


6.43180e5 


-6.90947 


1.35592e-7 


-5.03063e-4 


9.69105e-8 


0.00287 


46.83 


44.42 


340.89 


289.01 


85 


8(0,8) 


7.04940e5 


-7.32436 


5.07846e-9 


-2.49051e-4 


3.65309e-9 


0.00259 


56.26 


53.84 


440.57 


264.04 


60 


9(0,9) 


7.66200e5 


-7.72425 


1.79164e-10 


-1.24756e-4 


1.29718e-10 


0.00237 


66.47 


65.34 


584.40 


482.36 


82 


10(1,6) 


7.80170e5 


-7.09124 


6.40735e-5 


1.97283e-5 


4.94654e-5 


0.00238 


68.92 


71.52 


584.33 


529.65 


90 


11(3,2) 


8.37278e5 


6.82720 


-1.18614 


-2.25503e-6 


0.08218 


0.00235 


79.38 


77.75 


703.59 


581.81 


82 


12(2,6) 


9.66106e5 


-7.72737 


5.33342e-4 


-5.60824e-7 


4.52572e-4 


0.00192 


105.71 


105.10 


984.80 


782.46 


79 


13(0,14) 


10.67211e5 


9.55483 


2.51932e-14 


4.47644e-6 


1.77416e-14 


0.00167 


129.01 


124.64 


1101.62 


830.39 


75 


14(3,6) 


11.46563e5 


-8.32963 


0.00282 


1.87631e-8 


0.00265 


0.00161 


148.91 


143.28 


1202.54 


966.97 


80 


15(3,7) 


12.18265e5 


-8.62869 


2.67745e-4 


6.57974e-9 


2.52229e-4 


0.00150 


168.13 


162.58 


1654.50 


1244.89 


75 




7 8 9 10 11 12 13 14 
Acoustic Radiation Dominated Modes 



15 




0.2 0.3 0.4 0.5 
Orifice ratio (b/a) 



0.7 



FIG. 7. (color online) Comparison of experimental and com- 
puted quality factors 



FIG. 8. (color online) Effect of the central orifice size on 
the acoustic radiation damping in a few selected modes of 
vibration. The modes chosen here are (m, n) = (0,0), (0,1) 
and (1,1), taken purely as representative samples for circular 
and diametral modes. 



TABLE III. Effect of orifice size and modal parameters (m, n) 
on the acoustic damping of the resonator 









Qac 


for various 


b/a 




m 


n 


0.0 


0.1087 


0.3 


0.5 


0.7 








116.83 


119.01 


120.33 


140.11 


311.99 





1 


233.41 


236.42 


294.70 


235.13 


474.85 





2 


247.54 


260.02 


322.05 


320.74 


511.47 





3 


205.10 


206.13 


243.33 


302.45 


297.01 


1 


1 


198.43 


214.78 


249.30 


303.04 


739.40 


1 


2 


157.47 


161.15 


199.10 


354.23 


882.70 



resonator area results in much smaller reduction in the 
effective radiating area. 



IV. CONCLUSIONS: 

In this paper, we have presented a procedure for com- 
puting acoustic damping in various modes of vibration 
of an annular plate with fixed outer edge using exact 
mode shapes. A solid circular plate is treated as just a 
special case with the inner radius set to zero. The re- 
sults for a solid circular plate are compared with those 
of Lamb where he used approximate mode shapes and 
gave expressions for acoustic radiation losses in the first 
two modes. We have shown that the difference in the ex- 
act mode shape and the approximate mode shape results 
in considerable difference in the Q-factor, particularly in 
the second mode. We have also computed the Q-factor 
associated with acoustic radiation for various modes of 
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a micromechanical drum resonator and compared the re- 
sults with published experimental values where the ex- 
perimental results report the net Q in various modes. 
We have shown that the contribution of Q ac to the net Q 
is dominant for higher modes (above the 5 th mode of vi- 
bration) and accounts for almost 80% of the net Q. This 
result shows that for two dimensional micromechanical 
resonators, one should use those modes of vibration that 
emit least amount of acoustic energy to obtain highest 
possible Q. 
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APPENDIX A: SOLUTION PROCEDURE TO FIND 

w(r,6) 

The modal solution w(r, 9) of the eqn. ([3]) can be 
obtained by superimposing the solutions w\(r,9) and 
W2(t, 9) of the equations, 



V 2 (w l (r,9))±0 2 w l (r,( 



0, 



(Al) 



where, i = 1,2. Assuming the solution of the form 
Wi(r,9) — Xi{ r )-Tpi{9) based on the separable variable 
method and substituting it in eqn. (|A1[) . we obtain, 



[ d 2 Xi{r) 1 dxijr) 
dr 2 r dr 



Xi(r) 



1 d 2 A(9) 
d9 2 



= 0. 



Due to the symmetric nature of the problem^, periodic- 
ity in the azimuthal angle is introduced by equating the 
second term in the above expression to the separation 
constant —n 2 , we get, 

1 d 2 i> t {9) 2 



MO) d9 2 



r d 2 Xt(r) 1 dxijr) 
dr 2 r dr 



Xi(r) 



(A3) 



Under this condition, the modal solution can be written 
as w{r,9) = (xiM + X2( r ))' t P(9)- The equation (IA3|) is 
decomposed into two Bessel's equations^, in \i an d the 
simple harmonic equation^! in ip(9) as follow, 



d 2 Xi(r) , ldXi(r) 



dr 2 



dr 



± [Ft 



Xi{r) = 0, (A4) 



and 



d 2 ^(9) 
d9 2 



■n 2 i){9) = 0, 



(A5) 



where, n is the number of nodal diameters and it takes 
a value of only positive integer, i.e., n = 0, 1, 2, .... The 
solution of equation (|A5I) is ip(6) — C cos (nO T e), where 
A and e are arbitrary constants (for simplicity e = 0). 
The two Bessel's equations, 



d 2 Xi(r) 1 dxi(r) 
dr 2 r dr 



+ (P 2 -^)Xi(r) = 0, (A6) 



and 



d 2 X 2(r) , 1 dx2(r) _ ^ 2 + n 

r 



Xa(r) = 0, (A7) 



dr 2 r dr 
have the solutions of the form2£, 

Xi (r) = A rnn J n (/3r) + B mn Y n (j3r), (A8) 

and 

X2 (r) = C mn I n (fir) + D rnn K n {fir) , ( A9) 

where J n is the Bessel function of the first kind, Y n is 
the Bessel function of the second kind and A mn and B rnn 
are arbitrary constants in the solution of Xi( r )- In anc i 
K n are the modified Bessel's function of first and second 
kind, respectively, and C mn and D mn are arbitrary con- 
stants in the solution of X2(r). Consequently, the final 
form of the modal solution w(r,9) is found as^ 

w(r,6) = [A mn J n (f3r) + B mn Y n (J3r) + 
C m rMM + D mn K n (j3r)]i>(6) = x{r)i>{0). (A10) 



APPENDIX B: SYSTEM OF LINEAR AND 
HOMOGENOUS EQUATIONS FOR THE FOUR 
CONSTANTS A B C in n 

AND D 

in n 



A mn J n (fia)+B mn Y n (Pa)+C mn I n (pa)+D mn K n (0a) = 0, 

(Bl) 



- 1 / 



pa 



Jn(/3a) - J (n+ i)(/3o) 



+ 



B m 
D, 



0a 



13a 



Y n {f3a)-Y (n+1) {0a) 
l n (0a) +I (n+1) (j3a) 
K n {0a) - K {n+1) (0a) 



0, (B2) 



0b) + B mn F 2 (v,n,0b) 

- C mn F 3 (v, n, 0b) - D mn F A (y, n, 0b) = 0, (B3) 

0i (v,n,Pb) + B mn (j)2 (y, n,0b) 

- C mn (j)z(u, n, 0b) - D mn (j>i(v, n, 0b) = 0, (B4) 

where, the functions F±, F2, F3, F4, <p±, (f>%, <p3 and <f>4 axe 
defined by eqns. (|Cljl - (|C8p as mentioned on the ap- 
pendix [cl 
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APPENDIX C: EXPRESSIONS FOR FUNCTIONS 

Fx, F2, -F3, -F4, 4>i,<f>2, ct>3 AND 04 



F x (v,n,fib) = 
J n {fib)-(l-v) 



F 2 (p,n,fib) = 
Y n (fib)-{l-v) 



F 3 (v,n,fib) = 
I n {fib) + (l-v) 



F 4 (v,n,fib) = 
K n (fib) + 



(fibf 



(CI) 



{fib) 2 



(C3) 



-^T 1 K n (fib) + -K (n+1) (fib) 



§\(y,n,fib) = 



n 2 (l - v) 
((3b) 2 



nJ n (fib) - (fib)J {n+1) (fib) 
[{n-l)J n (fib)-{fib)J {n+l) (fib)] (C5) 



fa(v,n,Pb) = nY n (fib) - (fib)Y (n+1) (fib) 

, [(n - l)Y n ([3b) - (fib)Y {n+1) ((3b)] (C6) 



n 2 (l - v) 



(fibf 



(C2) cf> 3 {y,n,fib) = nl n (fib) + (fib)I (n+1) (fib) 



n 2 (l - v) 

(fib) 2 



[(n - l)I n (fib) + (fib)I (n+1) (fib)] (C7) 



Mv,n t pb) = nK n (fib) ~ (fib)K {n+1) {fib) 
n 2 (l - v) 



(fibf 



[(n - l)K n (fib) - (fib)K {n+1) {fib)] (C8) 



(C4) 



APPENDIX D: DETERMINANT OF THE COEFFICIENT 
MATRIX OF EQNS. [BT1fB4l 



MPa) Y n ([3a) I n ([3a) K n ([3a) 

jfcMfia)- J {n+1) (fia) ^Y n (fia)-Y {n+1) (fia) f-J n (fia) + I {n+1) (fia) ^K n (fia) - K (n+1) (f3a) 

Fx(v,n,pb) F 2 (v,n,fib) -F 3 (v,n,fib) -F 4 (v,n,fib) 

(f)i(p,n,pb) (f> 2 (v,n,/3b) -<fe(>, n, fib) ~4n(v,n, fib) 



APPENDIX E: DEFINITION OF f3 mn 

The frequency of free vibration of the plate in fluid 
can be related to its natural frequency in vacuum and is 
written asr^ 



fwet — 



I dry 



(El) 



where, 



is called added virtual mass incremen- 



tal(AVMI) factor^ as 



fir, 



Tf 
T P 



pj_a_ 

1 p s d 



(E2) 



fimn is defined as the ratio of reference kinetic energy Tp 
of the fluid, due to the structural vibration, and that of 
the resonator Tp. T mn is the non-dimensional added vir- 
tual mass incremental (NAVMI) factor. Note the above 
NAVMI factor is for plate in contact with fluid only on 



one side (Fig. [TJ. NAVMI factors obtained must be dou- 
bled for plates in contact with a fluid on both sides. 
NAVMI factors are found to be given by, 



r — 

i rn.'n. 



1 d p s 

fimn O, pf J Q 



H 2 mn ( V )dr,, (E3) 



where 



H mn (ri) — A mn HAmn(r]) + B mn HBmn(T]) 

+ C mn Hcmn(il) + D mn Hp> mn (ri), (E4) 

in which, 



HAmn(v) = 



1 



A 2 — T] 2 

b 

~ a(X 2 -T] 2 ) 
b 

a(A 2 — rj 2 ) 



[\Jn(ri)J n +i(X) - rjJ n+ i(r])J n (X)} 



AJ„ Jn+i ^A^ 

-rj,J n+ i (v-^j Jn f A- 



(E5) 
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Hsmniv) 



1 



A 2 — rj 
b 

a(X 2 — r/ 2 ) 
b 

a(X 2 - rj 2 ) 



2 [XJ n (v)Y n+1 (\) - nJ n+1 (r))Y n (X)] 



/ b\ ( b 

-■qJn+i I rj- J Y n I A- 



(E6) 



H. 



Cmn 



(v) 



1 



A 2 + r/ 2 
b 

a (A 2 + ?/ 2 ) 
6 

" a(A 2 + r ? 2 ) 



[AJ„(ry)J n+ i(A) + r/J n+ i(r/)J„(A)] 



6 








a 








( 









(E7) 



-1 



A 2 + 77 2 

6 



[AJ n (r ? )if„ +1 (A)-7 ? J n+1 (7 7 )^„(A)] 



a(A 2 +r? 2 ) 
a(A 2 +7? 2 ) 



/ 6\ / 6 

XJ n iv-J K n+ i I A- 

?yJn+l f 7 ?" ) (a^ 



(E8) 



and the above stated equations are used for annular res- 
onators. For solid disk resonator, we have, 

H mn (r]) = A mn H Amn(v) + C mn Hcmn(jl)i (E9) 



HAmniv) 



A 2 — rj' 



H, 



Cmn 



l 



2 [77J n _i(?7)J n (A) - AJ n _i(A) J n [ri)], 

(E10) 



r [T?J n+ i(r?)J„(A) + XJ n (r])I n+ i{X)}. 



A 2 + ?f 1 

(Ell) 

The velocity potential 'aV at the plate fluid interface 
is obtained by using the above H mn function, which has 
eliminatable singularity at rj = X (Here A = /3a) 17 ' 18 , 



/•OC 

(a,0) = -a / H mn (if)J n 
Jo 



{•qa)drj. 



(E12) 



The derivation of velocity potential assumes the fluid 
such as air to be incompressible, inviscid newtonian fluid. 
The plate is considered placed in an annular aperture of 
an infinite rigid wall. The fluid motion, considered only 
due to the plate's vibration, is assumed to be irrotational. 
The plate and fluid are maintained at ambient isothermal 
conditions. Hence such velocity potential must satisfy 
Laplace equation. By compatibility conditions velocity 
of the fluid in contact with the plate assumes the plate 
velocity. No slip flows occur at the fluid-structure in- 
terface. Such uncoupled fluid-structure interaction prob- 
lem involves evaluation of reference kinetic energy of the 
fluid at the plate-fluid interface by using the above ve- 
locity potential and hence this kinetic energy in-turn is 
used to evaluate /3 mn . The integral in eauation (IE3p was 
numerically evaluated in MATLAB using the adaptive 
Gauss-Kronrod quadrature. 
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